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Abstract 

The method of parameter variation for linear differential equations is extended to 
classes of second order nonlinear differential equations. This allows to reduce the 
latter to first order differential equations. Known classical equations such as the 
Bernoulli, Riccati and Abel equations are recovered in illustrated relevant examples. 
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1 Introduction 



Second order nonlinear differential equations are important for investigation 
of nonlinear phenomena in different fields of physics and mathematics. They 
are used to model a wide number of phenomena in plasma physics, solid state 
physics, optics, bio-hydrodynamics, chemical processes, nonlinear quantum 
mechanics, etc. See [T]|5] and references therein. Formal methods of analytical 
integration of nonlinear differential equations thus appear of great interest in 
the theory of differential equations. To cite a few, one can mention various 
powerful methods, which are most familiarly used, such as group symme- 
try methods [MI], tanh method [HlfTU], extended tanh method fTTlfT3] . sine- 
cosine method [T4lfT5] . Jacobi elliptic method [T6lfT7] . Backlund transforma- 
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tions p~8lTT9| . inverse scattering method [20], pseudo spectral method [21] and 
F-expansion method [22]|!M] . 



Unfortunately, only a few types of nonlinear differential equations can be ex- 
actly solved, what explains the permanent need to recur to novel tools of 
handling differential equations. In such a direction, the method of reducing 
the order of ordinary differential equations (ODEs) can be of inestimable use- 
fulness in most cases. Indeed, in general, the lowest the order of a differential 
equation is, the greatest is the possibility of finding suitable analytical methods 
for its resolution. Moreover, the accuracy of numerical treatment of a differen- 
tial equation decreases as the order of the equation increases. Therefore, when 
a differential equation cannot be directly analytically solved, it is meaningful 
to transform it into an equation of lower order. The probability of solving 
the reduced equation is greater than the probability of solving the original 
equation. In practice, it is generally difficult to perform such a reduction. 

This paper aims at proving that the usual method of parameter variation can 
be successfully used in some cases to reduce by one the order of nonlinear 
differential equations (NLDEs). In the sequel, we provide with some classes 
of second order nonlinear differential equations, transformable to first order 
ones using specific parameter variations. Concrete examples are exhibited as 
matter of illustration. Specifically, the classes of equations of the forms 



where m is a positive integer and a is an integrable function, are considered. 
Note that, if the function f(x,y,y') = g(x,y'), i.e. the second members of 
equations and (j2J) do not explicitly depend on the variable y, then we can 
perform the natural change of variables z(x) = y'(x) to lower the order of 
these equations. Besides, if the function a is constant only in equation (JTJ, 
(not necessarily constant in equation (j2j)), and if f(x,y,y') = g(y,y'), i.e. the 
equations ([1]) and (|2D are autonomous, then we can perform the change of 
variables w(y) = y'(x) to transform the above mentioned classes of equations 
into first order ones. 





and 




(2) 
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2 First class of reducible second order NLDEs 



Let us consider nonlinear second order differential equations of the following 
type: 

y" + a{x)y' = F(x) + y'G(y)e- 1 a ^ dx (3) 

whose the linear part, i.e. 

y" + a(x)y' = (4) 

yields the solution 

y > = Ce-f a{x)dx , (5) 

where C is an arbitrary constant. Suppose that C is a differentiate le function 
of both variables x and y expressed in the form 

C = H(x) + K(y). (6) 
Then, (E]) can be rewritten as 



y' 



[H{x) + K{y)]e-I a(x)dx (7) 



that we differentiate to obtain 



y" — g~ / a(x)doL 



H' + y'K 



a(x)y 



where 

H-^Jp and K= d JM, 

ax y ay 

Substituting Q and © into ©, we find 

H' x + y'^ = F(x)el a ^ dx + y'G(y). (9) 

Clearly, ((9]) will take place if 

H' x = F(x)el a{x)dx and K' y = G{y). (10) 

We therefore state the following result: 
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Proposition 2.1 The second order nonlinear differential equation can be 
reduced to the first order differential equation ([?]), where the functions H and 
K are solutions of the first order differential equations (fiOj). respectively. 

As a matter of clarity, let us consider the following example: 

Example 2.1 Let us consider the function G in ([3]) in the form 

G{)= dfb + bl y + b 2 y> + b 3 ys\ 
dy \ c + c x y J 

where b , b x , b 2 , b 3 are arbitrary constants and c , Ci are constants such that 
(c ,Ci) 7^ (0,0). Then equations (JTUJ) yield : 



H(x)= J F(x)el a{x)dx dx 



K(y)=A + bo + hy + hy2 + hy \ (12) 
c + ciy 

where A is an arbitrary constant of integration. Therefore, equation (j7j) is 
reduced to the well known Abel equation of second kind 



„/_«,- f &o + Cq(H(x) + A) + [ Cl (H(x) + A) + h] y + b 2 y 2 + b 3 y 3 

c + cxy 

(13) 

Finally, the equation (j3J) takes the form 

y" + a(x)y' = F(x) + y'±- ( b ° + b ^ + b ^ + b ^ ) e - (14 ) 

and is integrable if the constants bo, bx, b 2 , b 3 , Co, c\ and the functions a and 
F are chosen in such a way that the Abel equation (fl3|) be integrable. 



In particular, for Co = 1, c\ — and F — 0, equation ([TBI) leads to the 
separable Abel equation 

</ = (A + b + 6 l2/ + b 2 y 2 + 6 3 y 3 ) e-/^)* 8 (15) 



whose the implicit solution is given by 
f d V 



A + b + b iy + b 2 y 2 + b 3 y 3 



J e-I a{x)dx dx + B, (16) 



4 



where B is an arbitrary constant of integration. 

It is worth noticing that a second kind Abel equation of the form 

, hy Z + f2V 2 + fiy + fo , , n ( , 7 s 

V = ; , with f 3 0, (17) 

giy + go 

where f iy (i = 0, 1, 2, 3), and gj, (j = 0, 1), are arbitrary functions of x, can 
be transformed into a canonical form. Indeed, using the variable change 

L = t , y = l^\, (18) 

I giu J 

where t and m = u(t) are the new independent and dependent variables, 
respectively, equation (iTTj) becomes 

= / 3 n 3 + / 2 u 2 + /i« + fo. (19) 
Making use of the substitution 

u = v-A, (20) 

■J J 3 

equation flTTJT) can be put in the form 

= h 3 v 3 + h lV + ho. (21) 

Now, setting 

v = E(t)w, where £(t) = e/^ 1 ^', (22) 

brings this equation to the simpler form: 

w' t = h 3 w 3 + h , (23) 

which, in turn, can be reduced, with the help of the new independent variable 
h 3 (t)dt, (24) 



to the usual canonical form of Abel equation of the first kind 

w' s = w 3 (s) + k(s). (25) 
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The latter is integrable by various methods known in the literature. See [26] 
and [29] (and references therein) for a good compilation of techniques devel- 
oped to solve (1251) for particular expressions of k(s). 



3 Second class of reducible second order NLDEs 

In this section, we discuss the second relevant type of nonlinear second order 
differential equations compilable in the following general form: 

V" + a(y)(y'f = F(x)eS^ + y'G(y). (26) 

We first solve the left hand side part of this equation, namely 

y" + a(y)(yf = (27) 



to get 



y' 



Ce -Ja(y)dy^ ^ 



where C is an arbitrary constant. As in the previous section, we assume that 
C is a differentiable function of both x and y as follows 

C = H(x)+K(y). (29) 



Then, the equation (1281) takes the form 

y > = [H(x) + K{y)\ e~ f a(y)dy (30) 

that we differentiate to obtain 

'H' x + y'K' y ]-a(y)(y'f. (31) 



Substituting and ^ into ([26]), we find 

H' x + y'K' y = F(x) + y'G(y)eI a ^ d y. (32) 

Clearly, (!32|) will take place if 

H' x = F{x) and K' y = G(y)ef a{v)dy . (33) 
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Therefore, the following statement holds: 

Proposition 3.1 The second order nonlinear differential equation j[2h}) can be 

reduced to the first order differential equation / fffOj) . where the functions H and 
K are solutions of the two first order differential equations / fffffj) . respectively. 

For illustration, let us consider the following example. 

Example 3.1 Let a(y) = — i and G(y) = (3 y n , where n is a non zero positive 
integer and (3 is a constant. Then equation fliZoT) becomes 

y"-\y>y = F(x)y + {3y'y n . (34) 

y 

Equations (I3"3"j) yield 

H(x) = [ F(x)dx + A and K(y) = -y n } (35) 

where A is an arbitrary constant of integration. 

Therefore, equation ( |30l) becomes the Bernoulli equation [3]] 

= (#(*) + A) y + ^V 1 . (36) 



n 



The substitution w(x) = y l n transforms (1361) into the linear equation 

w ' x = - n (H(x) +A)w-p (37) 

whose the solution is 

. , -(3fe n I {H{x)+A)dx dx + B . . 

w[x) = 77- — - — - — , (38) 

V ' e nf (H(x)+A) dx V ' 

where B is an arbitrary constant of integration. 

4 Third class of reducible second order NLDEs 



The third group of second order nonlinear differential equations can be ex- 
pressed as 

{y') m y" + <x) {y') m+l = e^ m+2 ^ a ^ dx F (y, y' e I a ^ dx ) . (39) 
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Remark 4.1 If we set in Q3J F{x) = Ce~ 2 S a(x)dx , where C is a constant, 
then equation appears as a particular case of equations (1391) considered 
with m = 0. 

The linear part of equation (1391) 

(yT V" + a(x) (y') m+1 = (40) 

can be readily solved to give 

(y') m+1 = K m+1 e ~ (m+1) / a ( x ) dx ^ (41) 

where K is an arbitrary constant. Now, suppose that if is a differentiable 
function of the variable y. Then, the equation (141 p becomes 

(y') m+1 = K m+1 (y)e^ ( - m+1 ^ / a ^ dx , (42) 

Differentiate (142!) to obtain 

(2/)V n+2) = K m+1 (y)if;e'( m+2 )/ a W^ - a(x) {y') m+1 ■ (43) 

Substituting (PJ) and (PJ) into ((BSD, w e find 

K m+1 (y)^ = F(y,^)). (44) 

We therefore set the following result: 

Proposition 4.1 The second order nonlinear differential equation / fffgj) can 
be reduced to the first order differential equation ^J4\ ) using ^JW - Furthermore, 
if K(y) = $>(y, A), where A is an arbitrary constant, is the general solution 
of then the general solution of / fffPj) is given by 

I*£rrJ'- I "*' i * +B ' (45) 

where B is an arbitrary constant. 

The following examples are of particular significance. 

Example 4.1 Suppose in (1391) that F has one of the following forms 



i) F(u, v) = u m+1 f 



au + bv + c 

au + (3v + 7 
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or 

. .\ 7— t / \ m+ i , / cm + 6f + c \ . N 

u)F(u, v) = v m+1 f , 46 

yaw + pt> + 7/ 

where a, (3, 7, a, 6 and c are constants such that the function / is well de- 
fined. Then, using the method developed in [26] and [28], the corresponding 
first order differential equation (jUJ) can be transformed into a homogeneous 
equation which, in turn, can be reduced to a separable equation. As a matter 
of illustration, let us perform such a transformation in the case where F is of 
the form it). 

(1) If A = a/3 — ba ^ 0, the transformations 

67 — c(3 . . . . ca — cry , . 

2/ = ^ + A and ^ (3/) = + — ^— L (47) 

lead to the equation 

, „ / cm + bv \ , . 

< = / -75 • 48 

Dividing both the numerator and the denominator of the fraction on the 
right-hand side by u, we obtain the homogeneous equation 

for which the substitution u>(ii) = ^ gives the separable equation 

uw' u = J{w)-w. (50) 

(2) For A = and 6 7^ 0, the substitution t> (y) = ay + bK(y) + c engenders 
the separable equation 

f<9) For A = and (3 7^ 0, from the substitution t>(y) = ay + (3K(y) + 7, we 
deduce the separable equation 



'bv + c/3 — 67 



< = " + £ (52) 



Setting m = 0, a(x) = ± and F(it, v) = v (^) 2 + 2 (£) 
becomes 



the equation (1391 



.12 



y" + -y' = xy' [* ] +2 V -. (53) 
x \y I y 
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Perform in fl53l) the substitution (HZD which takes the form 



y = • (54) 

X 



By dI3|, the function K satisfies the homogeneous equation 



which gives the solution 

K{y) = (56) 
A-y 



where A is an arbitrary constant of integration. Hence, using the integral (145]) , 
we obtain, for equation fl53|) . the implicit solution 

\n(y) = \n(x) + B, (57) 

y 



where B is an integration constant. 

Suppose that F, defined in f l39l) . has the following general form 

F{u ' v) -EUk,(u)G v ( v y (58) 



where p, q are positive integers, k v , F v and are arbitrary functions such 
that F be well defined. 



Example 4.2 If F in (15811 is of the form 



F{u,v) = hx{u)R[v,JP{v)), (59) 



where R is a rational function of two variables and the function P under 
the radical is a polynomial of degree three or four, then the corresponding 
reduced equation (jSj) readily turns to be an elliptic integral which can be 
merely integrated by the method presented in [32] . 



Example 4.3 If F in (158[) has the form 

F{u,v) =v m+1 {h 1 {u)v + h 2 {u)v n ), n G IN \ {0,1}, (60) 
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then the corresponding reduced equation (PHI) leads to a Bernoulli equation. 
With m = 0, a(x) — 2/x and F(u, v ) = v (v + v 3 ) , the equation (I39I) becomes 



y" + -y'=(y') 2 + (xy') A . (61) 



.2' 



Perform in f[6"Tj) the substitution (1421) which takes the form 



x 2 



(62) 



By dB]) , the function satisfies the Bernoulli equation 

K' y = K + K 3 (63) 

yielding the solution 

^) = ±(Ae^-l) 4 , (64) 

where A is an arbitrary constant of integration. 

Hence, by using the integral fl58|) . the equation (|6T|) provides an implicit solu- 
tion 



arctan 



(vMe- 2 f - l) - VAe- 2 y - 1 = - + 5, (65) 



where B is an integration constant. 
Example 4.4 If F in fl58|) has the form 

F(u, v) = v m + l {h x {u) + h 2 (u) v + h 3 (u) v 2 ) } (66) 

then the corresponding reduced equation (jlij) leads to a Riccati equation. 
An important number of integrable Riccati equations is recorded in [2~TH3"U] . 
Recall that, if yo = yo(x) is a given particular solution of the Riccati equation 

y' = f(x)y 2 + g(x)y + h(x), (67) 



then, the general solution can be written as: 



y(x) = Vo{x) + *(x) 



C - J f(x)$(x)dx 



(6J 
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where 

$(x) = exp \ I [2f{x)y {x) + g{x)\ dx \ ; (69) 



C is an arbitrary constant. To the particular solution yo(x) there corresponds 
C = oo. 



Example 4.5 Let F in (joojl be of the form 

+1 + fe 2 («) v + h(u) v 2 + h A (u) v 3 

r{u,v)=v , , , , , , . UU 

fei(u) + fc 2 (u) u 



Then the corresponding reduced equation (jUJ) leads to an Abel equation. 

As a matter of fact, let's briefly present how to transform, into its canonical 
form, the Abel equations of the second kind: 

2/ = : , with / 2 ^ 0, (71) 

gi y + 9o 



where fa, i = 0, 1, 2, and pj, j = 0, 1, are arbitrary functions of x. 
Using the variable change 

y = (72) 

9i 

equation (17B takes the form 

zz' = hz 2 + hz + f . (73) 

Now, the substitution 

z = E(x) w where E{x) = e^ 2{x)dx (74) 

brings this equation to the simpler form : 

w w' x = hi w + ho, (75) 



which, in turn, can be reduced, by the introduction of the new independent 
variable 

hi(x)dx, (76) 
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into the canonical form of the second kind Abel equation 
w(s) w' s = w(s) + k(s). 



(77) 



A good compilation of integrable Abel equations of the form (j77J) can be found 



in 

Considering m — 0, a(x) = — - and F(u, v) = v + 2u, the equation (13^1) can 
be reduced to the linear equation 

y"-(l + x ) y'-2x 2 y = Q- (78) 



Perform in ( 1781) the substitution ( )42l) which takes now the form 

y' = xK{y). (79) 



By (H^l) . the function K verifies the Abel equation of second kind in its canon- 
ical form 

KK'=K + 2y. (80) 



Two particular solutions of equation (!80l are given by 

K l (y) = 2y and K 2 (y) = -y (81) 

from which can be deduced a general solution satisfying the algebraic equation 
(K(y)-2yf (K{y)+y)=A, (82) 



where A is an arbitrary constant. A real solution of the latter can be computed 
to yield 



K(y) = - y-8y 3 + AA + A^J-AAy 3 + A 2 

M v ^ 

+— ■ +y 

$-8 y 3 + 4 A + 4 V-4 Ay' s + A 2 
= $(y, A). 
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Hence, by using the integral (15*81) . we obtain an implicit solution of the equation 
781) as 



J [<%, AT^^J + B, (83) 
where B is a constant of integration. 

5 Fourth class of reducible second order NLDEs 



Let us now investigate second order nonlinear differential equations of the 
following form: 

(y') m y" + a{y) (y') m+2 = e" (m+1) J ' a{y)dy F (x, y'el a( - y)dy ) , (84) 

Remark 5.1 It is immediate to note that, if we set in (1261) G(y) = C, where 
C is a constant, then equation (126]) appears as a particular case of equation 
(1841) with m = 0. 



We first solve, as in the previous sections, the left hand side part of equation 
( 1841) . namely 

(yT y" + a(y) (y') m+2 = (85) 



from which follows 

(y') m+1 = K m+1 e~ {m+1) I a ^ dy , (86) 

where K is an arbitrary constant. Suppose that K is a differentiable function 
of the variable x. Then, equation (1861) becomes 

( y ') m+1 = K m+1 (x)e- {m+1) f a{y)dv . (87) 

Differentiate (1571) to obtain 

(y') m y " = K m (x)K' x e- {m+1) I a{y)dy - a(y) (y') m+2 ■ (88) 

Substituting (JEZD and flES) into §S) we find 

K m (x)K' x = F(x,K(x)). (89) 

We therefore arrive at the following result: 
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Proposition 5.1 The second order nonlinear differential equation can 
be reduced to the first order differential equation / fgpj) using ( fffTl ). Furthermore, 
if K(x) = $(x, A), where A is an arbitrary constant, is the general solution 
of / fffP|) . then the general solution of (5^] ) is given by 

J el a(y)dy dy = J $(x, A)dx + B, (90) 



where B is an arbitrary constant. 
Remark 5.2 

(1) The equations (13"9"j) and (I54"l) have been investigated by Jovan in [25] /or 
m = 0. Furthermore, this author has also examined some examples of 
integrable equations corresponding to particular cases where F is of the 
form (|3E!) with q = 1, fci(u) = = 1 and F v (v) = v a ", with a„ G M, 

i/ = 1, . . . , p. 

The equations (1551 and flM|) . considered with m ^ 0, are equivalent to 
the equations of the same type with m — 0. 

The examples given in the previous section can be also considered in the frame- 
work of the Proposition 15. 11 Namely, one can easily modify those examples to 
obtain equations of the form 0841) which can be integrated by quadrature in 
simpler situations or reduced at least to first order equations in more cumber- 
some cases. 

As a matter of illustration, let us consider the following examples in which the 
function F is of the form (1581) and m — 0. 



Example 5.1 Let a(y) = —1/y and F(u, v) = J (1 — v 2 ) (1 — k 2 v 2 ) where 
k G H*. Equation (1841) becomes 



y - - {y ) 
y 



(91) 



Perform in (19T|) the substitution (IH7I) which takes the form 

y' = K(x)y. (92) 

By ( |89|) . the function K satisfies the elliptic integral 
r dK 



'(1 - K 2 ) (1 - A; 2 X 2 ) 



x + A, (93) 
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where A is an arbitrary constant of integration. Equation (193]) is readily solved 
to give 

K{x) = sn(x + A, k), (94) 

where sn is the first Jacobi elliptic function. Hence, by using the integral (I9"U|) . 
the solution of the equation fl9~Tl) is obtained in the form 

y(x) = B [dn(x + A, k) + k cn(x + A, k)]~* , (95) 

where B is an integration constant; cn and dn stand for the second and third 
Jacobi elliptic functions, respectively. 



Example 5.2 Let a(y) = and F(u, v ) = + av 2 , a, b G 1R*. Equation 
becomes 



y" - - (y'T = y 
y 



b , iy 

-tt + a - 



x- 



y 



(96) 



Perform in (196]) the substitution (IH7I) which takes the form 
y' = K(x)y. 



(97) 



By fl89|) . the function K verifies the Riccati equation 



K' x = — + aK 2 



(9£ 



which can be solved to yield 



2aX+l \ - 1 



(99) 



where A is an arbitrary constant of integration and A is a solution of the 
quadratic equation a A 2 + A + b = 0. Hence, by using the integral (}9~D1) . we 
obtain, for equation (1961) . the solution 



(aie 2aAlnx + 2aAA + A)' 



(100) 
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where B is a constant of integration. 



6 Concluding remarks 

We have investigated, in this paper, four different classes of second order non- 
linear differential equations which have been reduced to first order ones, using 
suitable parameter variations. Fortunately, the resulting first order differential 
equations are, in most cases, transformable to well known integrable or solv- 
able classical differential equations whose the solutions can be worked out by 
various methods disseminated in the standard specialized text books. Finally, 
it appears possible to extend the parameter variation methods developed in 
this work to classes of higher order nonlinear differential equations with a 
view to their order reduction. Such an investigation will be in the core of 
forthcoming work. 
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